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ABSTRACT 

In this Paper Hybrid techniques are the integration of 
various soft computing techniques that synergizes the 
strength and weakness of an individual technique with 
the other therefore the overall performance of such 
techniques working in unison is increased. A 
deterministic inventory of Auto industry model has 
been developed for deteriorating items of Auto 
industry having a ramp type demand with the effects 
of inflation with two-storage of Auto industry 
facilities using Particle Swarm Optimization. Here, 
we assumed that the inventory of Auto industry 
holding cost in RW is higher than those in OW. 
Shortages in inventory of Auto industry are allowed 
and partially backlogged and it is assumed that the 
inventory of Auto industry deteriorates over time at a 
variable deterioration rate using Particle Swarm 
Optimization. The effect of inflation of Auto industry 
has also been considered for various costs associated 
with the inventory of Auto industry system. Cost 
minimization technique is using Particle Swarm 
Optimization to get the expressions for total cost and 
other parameters. 

Keywords: Two-warehouses, deterministic inventory, 
deteriorating items and Particle Swarm Optimization 

1. INTRODUCTION 

Many researchers extended the EOQ model to time- 
varying demand patterns. Some researchers discussed 
of inventory of Auto industry models with linear trend 
in demand. The main limitations in linear-time 
varying demand rate is that it implies a uniform 
change in the demand rate per unit time. This rarely 
happens in the case of any commodity in the market. 


In recent years, some models have been developed 
with a demand rate that changes exponentially with 
time. For seasonal products like clothes, Air 
conditions etc. at the end of their seasons the demand 
of these items is observed to be exponentially 
decreasing for some initial period. Afterwards, the 
demand for the products becomes steady rather than 
decreasing exponentially. It is believed that such type 
of demand is quite realistic. Such type situation can 
be represented by ramp type demand rate. 

An important issue in the inventory theory is related 
to how to deal with the unfulfilled demands which 
occur during shortages or stock outs. In most of the 
developed models researchers assumed that the 
shortages are either completely backlogged or 
completely lost. The first case, known as backordered 
or backlogging case, represent a situation where the 
unfulfilled demand is completely back ordered. In the 
second case, also known as lost sale case, we assume 
that the unfulfilled demand is completely lost. 

Furthermore, when the shortages occur, some 
customers are willing to wait for backorder and others 
would turn to buy from other sellers. In many cases 
customers are conditioned to a shipping delay and 
may be willing to wait for a short time in order to get 
their first choice. For instance, for fashionable 
commodities and high-tech products with short 
product life cycle, the willingness of a customer to 
wait for backlogging is diminishing with the length of 
the waiting time. Thus the length of the waiting time 
for the next replenishment would determine whether 
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the backlogging would be accepted or not. In many 
real life situations, during a shortage period, the 
longer the waiting time is, the smaller is the 
backlogging rate would be. Therefore, for realistic 
business situations the backlogging rate should be 
variable and dependent on the waiting time for the 
next replenishment. Many researchers have modified 
inventory policies by considering the “time 
proportional partial backlogging rate”. 

2. RELATED WORKS 

Zangwill (1966) developed a production multi period 
production scheduling model with backlogging. 
Inventory models with a mixture of backorders and 
lost sales were formulated by Montgomery et al. 
(1973). Economic production lot size model for 
deteriorating items with partial backordering was 
suggested by Wee (1993). A comparison of two 
replenishment strategies for the lost sales inventory 
model was presented by Donselaar et al. (1996). 
Time-limited free back-orders EOQ model was 
formulated by Abbound and Sfairy (1997).Chang and 
Dye (1999) developed an inventory model in which 
the proportion of customers who would like to accept 
backlogging is the reciprocal of a linear function of 
the waiting time. Papachristos and Skouri (2000) 
established a partially backlogged inventory model in 
which the backlogging rate decreases exponentially as 
the waiting time increases. An EOQ inventory model 
for items with Weibull distribution deterioration rate 
and ramp type demand was formulated by Wu (2001). 
Shortages in inventory were allowed and assumed to 
be partially backlogged. Teng et al. (2002) presented 
an optimal replenishment policy for deteriorating 
items with time-varying demand and partial 
backlogging. An EOQ model for deteriorating items 
with time-varying demand and partial backlogging 
was suggested by Teng et al. (2003). Chu and Chung 
(2004) discussed the sensitivity of the inventory 
model with partial backorders. An EOQ model with 
time varying deterioration and linear time varying 
demand over finite time horizon was proposed by 
Ghosh and Chaudhuri (2005). Shortages in inventory 
were allowed and partially backlogged with waiting 
time dependent backlogging rate. Optimal ordering 
policy for deteriorating items with partial backlogging 
was formulated by Ouyang et al. (2006) when delay in 
payment was permissible. Dye (2007) proposed joint 
pricing and ordering policy for deteriorating 
inventory. Shortages in inventory were allowed and 
partial backlogged. An inventory lot-size model for 
deteriorating items with partial backlogging was 


formulated by Chem et al. (2008). Authors have taken 
time value of money in to consideration. The demand 
was assumed to fluctuating function of time and the 
backlogging rate of unsatisfied demand was a 
decreasing function of the waiting time. 

Taleizadeh AA, et.al. (2008). developed Optimizing 
Multi-Product Multi-Constraint Inventory Control 
Systems with Stochastic Replenishment and an EOQ 
Joint Replenishment Policy to Supply Expensive 
Imported Raw Materials with Payment in Advance, 
multi products, multi constraints newsboy problem 
with batch order and discount. 

Taleizadeh AA, et.al (2009) extended A Particle 
Swarm Optimization Approach for Constraint Joint 
Single Buyer-Single Vendor Inventory Problem with 
Changeable Lead-Time and (r,Q) Policy in Supply 
Chain, developing a Multiproducts Newsboy Problem 
with a Hybrid Method of Fuzzy Simulation and 
Simulated Annealing, a Hybrid Method of Pareto, 
TOPSIS and Genetic Algorithm to Optimize Multi- 
Product Multi-Constraint Inventory Control Systems 
with Random Fuzzy Replenishments, optimizing 
multi product multi constraints Bi-objective newsboy 
Problem with discount by Hybrid Method of Goal 
Programming and Genetic Algorithm, a Hybrid 
Method to Optimize Bi-Objective Single Period 
Newsboy Problem with Fuzzy Cost and Incremental 
Discount, multi-product multi-constraint inventory 
control systems with stochastic replenishment and 
discount under Fuzzy Purchasing Price and Holding 
Costs, multi-product multi-constraint inventory 
control systems with stochastic period length and total 
discount under fuzzy purchasing price and holding 
costs, optimizing multi-product multi-constraint bi¬ 
objective newsboy problem with discount by hybrid 
method of goal programming and genetic algorithm, a 
harmony search algorithm to optimize chance 
constraint stochastic inventory systems with dynamic 
demand, joint replenish-up-to multi products EOQ 
with fuzzy rough demand. N. Jawahar and A.N. Balaji 
(2009) developed A genetic algorithm for the two- 
stage supply chain distribution problem associated 
with a fixed charge. Haixin Zhang et.al. (2013) 
suggested A modified multi-criterion optimization 
genetic algorithm for order distribution in 
collaborative supply chain. S.P. Nachiappan and N. 
Jawahar (2007) developed A genetic algorithm for 
optimal operating parameters of VMI system in a two- 
echelon supply chain. Z.H. Che and C.J. Chiang 
(2010) suggested a modified Pareto genetic algorithm 
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for multi-objective build-to-order supply chain 
planning with product assembly. Keyvan Sarrafha 
et.al. (2015) developed A bi-objective integrated 
procurement, production, and distribution problem of 
a multi-echelon supply chain network design: A new 
tuned MOEA. Ata Allah Taleizadeh et.al. (2011) 
extended Multiple-buyer multiple-vendor multi¬ 
product multi-constraint supply chain problem with 
stochastic demand and variable lead-time: A harmony 
search algorithm. Wei-Chang Yeh and Mei-Chi 
Chuang (2011) suggested Using multi-objective 
genetic algorithm for partner selection in green supply 
chain problems. S.A. Torabi, et.al. (2006) extended A 
hybrid genetic algorithm for the finite horizon 
economic lot and delivery scheduling in supply 
chains. Pavan Kumar Naraharisetti et.al. (2007) 
suggested Optimal supply chain redesign using 
genetic algorithm. Moon-Chan Kim et.al. (2008) 
developed Forward-backward analysis of RFID- 
enabled supply chain using fuzzy cognitive map and 
genetic algorithm. Partha Guchhait et.al. (2013) 
extended A production inventory model with fuzzy 
production and demand using fuzzy differential 
equation: An interval compared genetic algorithm 
approach. U.K. Bera et.al. (2012) extended Inventory 
model with fuzzy lead-time and dynamic demand 
over finite time horizon using a multi-objective 
genetic algorithm. Ata Allah Taleizadeh et. al. (2009) 
suggested A hybrid method of Pareto, TOPSIS and 
genetic algorithm to optimize multi-product multi¬ 
constraint inventory control systems with random 
fuzzy replenishments. Ali Roozbeh Nia et.al. (2014) 
extended A fuzzy vendor managed inventory of multi¬ 
item economic order quantity model under shortage: 
An ant colony optimization algorithm. Javad Sadeghi 
et. al. (2014) extended optimizing a hybrid vendor- 
managed inventory and transportation problem with 
fuzzy demand: An improved particle swarm 
optimization algorithm. Manas Kumar Maiti (2011) 
developed A fuzzy genetic algorithm with varying 
population size to solve an inventory model with 
credit-link ed promotional demand in an imprecise 
planning horizon. H. Altay Guvenir and Erdal Erel 
(1998) suggested Multicriteria inventory classification 
using a genetic algorithm. Ali Diabat (2014) extended 
Hybrid algorithm for a vendor managed inventory 
system in a two-echelon supply chain. K.L. Mak et. 
al. (1999) developed optimal inventory control of 
lumpy demand items using genetic algorithms. A.K. 
Maiti et. al. (2006) extended An application of real- 
coded genetic algorithm (RCGA) for mixed integer 
non-linear programming in two-storage multi-item 


inventory model with discount policy. Kuo-Ping Lin 
et.al. (2010) developed A simulation of vendor 
managed inventory dynamics using fuzzy arithmetic 
operations with genetic algorithms. Sung-Pil Hong 
and Yong-Hyuk Kim (2009) suggested A genetic 
algorithm for joint replenishment based on the exact 
inventory cost. Manas Kumar Maiti and Manoranjan 
Maiti (2007) extended Two-storage inventory model 
with lot-size dependent fuzzy lead-time under 
possibility constraints via genetic algorithm. Arindam 
Roy et. al. (2009) developed a production inventory 
model with stock dependent demand incorporating 
learning and inflationary effect in a random planning 
horizon: A fuzzy genetic algorithm with varying 
population size approach. S.M. Disney et. al. (2000) 
suggested Genetic algorithm optimisation of a class of 
inventory control systems. 

3. Assumptions and Notations 

In developing the mathematical model of the 
inventory system the following assumptions are being 
made: 

> A single item is considered over a prescribed 
period T units of time. 

> The demand rate of Auto industry D(t) at time t is 
deterministic and taken as a ramp type function of 

time i.e. D(t) = T)e~ r ^~ (t ~' l/>H(t ~ v) ^ 1 77 > 0, y > 0 

> The replenishment rate is infinite and lead-time is 
zero. 

> When the demand for goods is more than the 
supply. 

> The variable rate of deterioration in both 
warehouse is taken as P (t) = pt. Where 0< P « 1 
and only applied to on hand inventory. 

> No replacement or repair of deteriorated items is 
made during a given cycle. 

> The owned warehouse of Auto industry (OW) has 
a fixed capacity of W units; the rented warehouse 
of Auto industry (RW) has unlimited capacity. 

> The goods of OW are consumed only after 
consuming the goods kept in RW. 

In addition, the following notations are used 
throughout this paper: 

low(t) The inventory of Auto industry level in OW at 
any time t. 

I rw (t) The inventory of Auto industry level in RW at 
any time t. 

Cw The capacity of the own warehouse of Auto 
industry. 

Q The ordering quantity of Auto industry per 
cycle. 
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T Planning horizon. 

6 Inflation rate of Auto industry. 

Ho The holding cost of Auto industry per unit per Ci 
unit time in OW. 

Hr The holding cost of Auto industry per unit per A 
unit time in RW. where Ho < Hr 

Dc The deterioration cost of Auto industry per 
unit. 

4. Formulation and solution of the model 

The inventory levels of Auto industry at OW of Auto industry are governed by the following differential 

equations: 


The shortage cost of Auto industry per unit per 
unit time. 

The opportunity cost of Auto industry due to 
lost sales. 

The replenishment cost of Auto industry per 
order. 


^_ow(0_ 

dt 

dt 

And 

dJ 0W (t) 

dt 


= -P(t)i(t), 

+ P(t)I(t) = -ne- rv , 


-at 

= -tje ' Y e 


0 < t < iff 
y/ <t<t x 

tf < t < T 


( 1 ) 

( 2 ) 

(3) 


with the boundary conditions, 
If)w(O) =Cw and I(ti) = 0 


The solutions of equations (1), (2) and (3) are given by 
Iow0)-Cw e ^ ^, 0 <t<y/ 

i„J0= L v <, < (| 


And I ow d) = !le-"{e- a '-e - a «}. 
respectively. 


tj <t<T 


(4) 

(5) 

( 6 ) 

(7) 


The inventory level of Auto industry at RW of Auto industry is governed by the following differential 
equations: 


dIrw(t) - + j3(t)I(t) = -r 1 e- rt 


dt 


0 <t<y/ 


( 8 ) 


With the boundary condition I rw (0)=0, the solution of the equation (8) is 

P ^ ^ /? ' 2 

2 V > ~6 


Irw( t ) = v\(¥-t)~[Y 2 + )\e pt2 / 2 , l//<t<t x 


(9) 


Due to continuity of I 0 (t) at point t = p, it follows from equations (5) and (6), one has 


C w e~^ ' 2 = 


7je 


' ry , \( h - w ) + P ( ^ ^ 3 ) \ e ~^ 2 / 2 


C w -rj e w \ (ty -y/) + p 


($ V; 


( 10 ) 
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The total average cost consists of following elements: 

I. Ordering cost of Auto industry per cycle = A 

II. Holding cost of Auto industry per cycle (Hi) in OW 

¥ h 

O J Io w (t)e~ 9t dt + | I 0w (t)e~ e ^ +t ^dt 

0 y, 


H , =H, 


( 11 ) 


( 12 ) 


III. Holding cost of Auto industry per cycle (Chr) in RW 


Ht=H 


R 


¥ 


J I nv (t)e~ dt dt 


0 


(13) 


IV. Cost of deteriorated units of Auto industry per cycle (C D ) 

¥ ¥ h 

'C J Ptl rw (t)e~ 6t dt+ ^ptI Qw (t)e~ 6t dt+^ptI Qw (t)e~ 9(t+y/) dt 
0 0 y/ 


= D, 


(14) 


V. Shortage cost of Auto industry per cycle (Cs) 
C \- I 0w( t ) e ~ 6(h+t)dt 

A 


= S , 


(15) 


VI. Opportunity cost of Auto industry due to lost sales per cycle (C 0 ) 


1 

C x ^ri(\-e at )e w e~ d(h+t> dt 


(16) 


Therefore, the total average cost per unit time of our model is obtained as follows 
K(t,,T) = ^[Ordering cost of Auto industry + Holding cost of Auto industry in OW+ Holding cost of Auto 
industry in RW+ Deterioration cost of Auto industry + Shortage cost of Auto industry + Opportunity cost of 
Auto industry ] = ^ (say) (17) 
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5. Particle Swarm Optimization (PSO) 

Particle Swarm Optimization (PSO) introduced by Kennedy and Eberhart in 1995 is a population based 
evolutionary computation technique. It has been developed by simulating bird flocking fish schooling or 
sociological behavior of a group of people artificially. Here the population of solution is called swarm which is 
composed of a number of agents known as particles. Each particle is treated as point in d-dimensional search 
space which modifies its position according to its own flying experience and that of other particles present in 
the swarm. The algorithm starts with a population (swarm) of random solutions (particles). Each particle is 
assigned a random velocity and allowed to move in the problem space. The particles have memory and each of 
them keeps track of its previous (local) best position. 


U b 

Y b 

b 

q 

m 

ai, a 2 
e and 

nb 

pbest 

nb 

gbest 


Particle’s velocity 
Particle’s position 
number of elements in a particle, 

Inertia weight of the particle, 
generation number, 
acceleration constants, 
random value between 0 and 1 
Local best position of the particle, 

Global best position of particle in the swarm 


U b =qx [/*_! + % x eandi x [U b best - Y b _ r \ + a 2 x eand 2 x [U b best - Y b _ ± ] 

yb — V b 4 - II b 
1 n 1 n-1 ^ u n 
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6. Numerical Illustration 

To illustrate the model numerically the following parameter values are considered. 


r| = 50 units, 
0 = 0.05 unit, 
y = 0.2 unit, 

P =0.002 unit, 
Y = 0.2 year, 
a = 0.1 unit, 


A= Rs. 100 per order, 

Ho = Rs. 3.0 per unit per year, 
Hw = Rs. 10.0 per unit, 

Sc = Rs. 12.0 per unit per year, 
Ci = Rs. 4.0 per unit, 

T = 1 year, 


Then for the minimization of total average cost and with help of software, the optimal policy can be obtained 
such as: 

ti = 0.799224 year, S = 38.597235 units and K = Rs.158.115354 per year. 

Therefore, numerical solution of these equations is obtained by using the software MATLAB 7.0.1, 

Particle Swarm Optimization (PSO) 

population= 90, generations=600, cognitive learning factor=4, cooperative factor=4, social learning 
factor=0.10, inertial constant=0.10 and number of neighbors=10. 

Since the model is integer in nature, reaching an analytical solution (if any) to the problem is difficult (Gen and 
Cheng, 1997). So we need to use Meta heuristic algorithms. To solve the models under meta-heuristic 
approach, four hybrid intelligent algorithms of harmony search (Taleizadeh et al., 2009), simulated annealing 
(Taleizadeh, 2008, 2009), genetic algorithm Taleizadeh (2008, 2009), and particle swarm optimization 
(Taleizadeh et al., 2009) are used. A comparison of the results in Tables 1 and 2 for increasing and decreasing 
demand show that Hybrid Algorithms method performs the best. 


Table: - 1 Best results of purpose function by different algorithms for Maximum demand 


Hybrid Algorithms 

Product’s Maximum inventory 

Level 

Maximum 


1 ■■ 2 ■■ 3 

4 

5 

profit (Rs.) 

Particle Swarm Optimization 

720 803 370.4 

522.5 

782 

3012.15 


Table: - 2 Best results of purpose function by different algorithms for Minimum demand. 


Hybrid Algorithms 

Product’s Maximum inventory! 

Level 

Maximum 


1 

2 

3 

4 

5 

profit (Rs.) 

Particle Swarm Optimization 

520 

603 

170.4 

322.5 

582 

2012.15 


7. Conclusion 

The main pillars of soft computing are Particle Swarm 
Optimization all of these have been successfully 
applied in isolation to solve practical problems in 
various fields where conventional techniques have 
been found inadequate. However these techniques are 
complementary to each other and they may work 
synergistically combining the strengths of more than 
one of these techniques as the situation may 
occasionally demand. Hybrid systems are systems 
where several techniques are combined to solve a 
problem. Needless to say that such amalgamation 
must be used only if it returns better results than any 
of these techniques in isolation. Based on how or 
more systems are combined hybrid systems have been 


classified into three broad categories viz sequential, 
auxiliary and embedded hybrid systems. 

This study incorporates some realistic features that are 
likely to be associated with the inventory of any 
material. Decay (deterioration) overtime for any 
material product and occurrence of shortages in 
inventory are natural phenomenon in real situations. 
Inventory of Auto industry shortages are allowed in 
the model using Particle Swarm Optimization. In 
many cases customers are conditioned to a shipping 
delay, and may be willing to wait for a short time in 
order to get their first choice. Therefore, this concept 
is also taken in this model. From the numerical 
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illustration of the model, it is observed that the period 
in which inventory holds of Auto industry increases 
with the increment in backlogging and ramp 
parameters while inventory period decreases with the 
increment in deterioration and inflation parameters of 
Auto industry using Particle Swarm Optimization. 
Initial inventory level of Auto industry decreases with 
the increment in deterioration, inflation and ramp 
parameters while inventory level increases of Auto 
industry with the increment in backlogging parameter 
using Particle Swarm Optimization. The total average 
cost of the system goes on increasing with the 
increment in the backlogging and deterioration 
parameters while it decreases with the increment in 
inflation and ramp parameters using Particle Swarm 
Optimization. The proposed model can be further 
extended in several ways. For example, we could 
extend this deterministic model in to stochastic model. 
Also, we could extend the model to incorporate some 
more realistic features, such as quantity discount or 
the unit purchase cost, the inventory holding cost and 
others can also taken fluctuating with time using 
Particle Swarm Optimization. 
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